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Two-component repulsive Fermi gases with population imbalance in elongated 
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We study the two-component repulsive Fermi gas with imbalanced populations in one dimension. 
Starting from the Bethe Ansatz solution we calculate analytically the phase diagram for the homo- 
geneous system. We show that three phases appear: the balanced phase, the fully polarised phase 
and the partially polarised phase. By means of the local density approximation and the equation 
of state for the homogeneous system we calculate the density profile for the harmonically confined 
case. We show that a two-shell structure appears: at the center of the cloud we find the partially 
polarised phase and at the edges the fully polarised one. The radii of the inner and outer shells 
are calculated for different values of the polarisation and the coupling strength. We calculate the 
dependence of the magnetisation on the polarisation for different values of the coupling strength 
and we show that the susceptibility is always finite. 

PACS numbers: 03.75.Ss, 05.30.Fk, 71.10.Pm, 03.75.Mn 



I. INTRODUCTION 

The experimental achievement of a quantum degen- 
erate ultracold Fermi gas [H, 0, Q makes it possible to 
explore experimentally many new models, allowing us to 
study the rich physics of fermions At the same time, 
recent experiments have shown the possibility to study 
ultra-cold atomic gases confined in very elongated traps 
[H> El 0> II] ) stimulating theoretical interest in the physics 
of low-dimensional systems. The one-dimensional regime 
is reached through tightly confining the radial motion of 
atoms to zero point oscillations [1 31 ] ; this can be done ex- 
perimentally by loading the atoms in a sufficiently strong 
optical lattice [8( or by means of magnetic wires in an 
atom chip 

Present facilities also allow to create spinor gases by 
confining in the same trapping potential different atoms 
or several different internal states of the same atom, al- 
lowing us to study the physics of multicomponent atomic 
systems @, 

In addition, the interatomic interactions can be con- 
trolled via a Feshbach resonance [10| , making it possible 
to explore different interaction regimes. 

Combining the one-dimensional geometry with the 
tunable interactions and the possibility of mixing differ- 
ent species of atoms allows us to explore different features 
such as the Tonks-Girardeau gas 11, 
confinement induced molecules \lc 
of the BEC-BCS crossover! 
one-dimensional systems [16 1. 



121, the formation of 



l4 |15[ , the physics 



1181 ] . or the dynamics of 
Currently, uniform one- 



dimensional systems are relatively well understood and 
described. However, the presence of the shallow trapping 
potential in the long direction of the ultra-cold atomic 
cloud introduces a source of inhomogeneity in the sys- 
tem, bringing to light new physics which is particularly 



important in the case of multicomponent atomic systems. 

Polarised Fermi gase s have been the focus of interest 
in the recent years 17, 18j]. In this article we present a 
study of the exact ground state at zero temperature of a 
polarised one-dimensional Fermi gas with repulsive inter- 
component interactions. We show that for the harmon- 
ically confined case the atoms are distributed in a two- 
shell structure: a partially polarised phase is localised in 
the inner shell and a fully polarised phase sits at the edges 
of the trap. The radii of the inner and outer shell are cal- 
culated as a function of the polarisation. We calculate 
the dependence of the magnetisation on the polarisation, 
showing that the system always has a finite susceptibility. 
We begin by reviewing the Bethe Ansatz solution for the 
Hamiltonian of the mixture of a two-component Fermi 
gas with a delta-function repulsive interaction; then we 
use this solution to calculate the phase diagram of the 
system in the homogeneous case. Adding a harmonic 
confinement and using local density approximation, we 
calculate the properties of the system in the trapped con- 
figuration. 



II. MODEL 

Let us consider an atomic cloud confined in an 
anisotropic harmonic potential 
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(i) 



where m is the mass of an atom and r — y ' x 2 + y 2 . 
When the ratio between the two trapping frequencies sat- 
isfies lu z /uj± <C 1, only the lowest transverse mode is pop- 
ulated; then the low-energy scattering properties of such 
a system can be modeled by an effective contact inter- 
action in one dimension U{z) — gn)8{z) with parameter 
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gx D = 2fi 2 /m a 1D = 2 h 2 a 3D /[ma 2 L (1 - Ca 3D /a±)} [lj], 
where a± = y / n/muT[ is the transverse oscillation length, 
am is the effective ID scattering length, C = 1.0326... is 
a numerical parameter and a^u is the three-dimensional 
scattering length. By using magnetic-field-induced Fes- 
hbach resonances, the interaction can be tuned from a 
repulsive effective interaction (giD > 0) to an attractive 
one (g 1D < 0). 

We consider a one-dimensional two-component system 
of spin- 1/2 fermions interacting via a delta- function po- 
tential. The Hamiltonian of the system is 



H : 




Zj)\+V a «{z), (2) 



where c is related to the ID effective interaction coef- 
ficient by c = mgio/Ti 2 = 2 /am an< 4 Vcxt is the con- 
fining potential. Due to the Pauli exclusion principle, 
the delta-function interaction only occurs between parti- 
cles with opposite spin. The total number of particles is 
N = N-\ + Ny, without loss of generality we can assume 
N\ < N/2. We will call the particles with up (down) spin 
the majority (minority) component. In one dimension 
the interactions are measured through a dimensionlcss 
coupling constant 7 = c/n, which can be interpreted as 
the ratio between the mean field interacting energy and 
the characteristic kinetic energy. 7 < 1 corresponds to 
the weakly interacting regime while 7 ^> 1 stands for the 
strong interacting one. Hence, in one dimension, the gas 
is more strongly interacting the lower the density. 

In the homogeneous case, without the harmon ic, c on- 
finement, this system was exactly solved by Yang [20j for 
the case of repulsive interaction using the Bethe Ansatz: 
an eigenstate of the Hamiltonian Eq. @ is characterised 
by a set of N density quantum numbers Ij (j = 1, . . . , N), 
and iVj spin quantum numbers J a (a = 1, . . . , N^); which 
define a set of TV quasi-momenta pj and N± spin rapidi- 
ties Aq,. In the thermodynamic limit, they satisfy the set 
of Bethe Ansatz coupled equations: 
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(3) 



where all the parameters have been rescaled by c: A = 
A/c, p = p/c, B — B/c, Q = Q/c. pip) should be under- 
stood as the density of p (density of particles) and cr(A) 
as the density of A (density of down-spin particles). The 
cut-offs B and Q are determined by fixing the density of 
down-spin particles (njj and the total density of particles 
(n = nj, + Tij) in the system: 



n f +Q n, f +B 

- = / p(p)dp, ^ = / ai\)d\ 

C J-O c J—B 



(4) 



The weak coupling regime corresponds to n/c > 1 and 
the strong coupling regime to n/c < 1. The magnetisa- 
tion is defined as s = n j — n j , and the expression for the 
energy reads 

he 3 f Q 

P 2 pip)dp. (5) 
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III. PHASE DIAGRAM 

We will now calculate the phase diagram for the sys- 
tem. In order to do so, we fix the total density of particles 
n and change the balance between the numbers of up and 
down spins, that is to say the magnetisation s. When 
s = 0, the system is at balance, since the number of up 
spins equals the number of down spins, rif = Hi. Since 
the excitation spectrum is gapless this phase only exists 
in the situation with zero magnetic field. For s = n the 
system is fully polarised and n_[ = 0. For any intermedi- 
ate value of the magnetisation < s < n the system is an 
imbalanced mixture of up and down spins, i.e. partially 
polarized. There is no evidence of a finite momenta in- 
stability, such as in the FFLO state , in the partially 
polarised phase: a calculation of the critical exponents 
for the pair correlation function shows that this correla- 
tion is never the leading one for the case with repulsion 

21 1 , contrary to what happens in the case of attraction 

221. 



The magnetic field and the chemical potential, 
dE(n,s)/L dEjn, s)/L 

k= ds ' " = On ' (6) 
can be calculated from Eq. ((3]), (f4]) and ((5|): the energy 
depends directly on the cut-offs Q and B, so one has to 
take the derivative of (|5|) with respect to Q and B and 
then use (U]) to find the derivatives of Q and B with 
respect to n and s. The magnetic field and chemical 
potential can also be related to the chemical potential of 
the two species of fermions as h — MT ~ MJ - and /1 = A ' T ^" A ' 1 . 

Setting s = n, or equivalently = in Eq. ([S]), one 
can calculate the boundary between the mixed imbal- 
anced phase, where < s < n, and the fully polarised 
phase corresponding to s — n. The saturation magnetic 
field and the corresponding chemical potential for which 
the system becomes fully polarised are 



h s = 2dB 

(i a = -h 
where e# = 



— arctan(2Qo) ~ tt~ 

7T " 2lT 

2e B Ql 
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h c : /4m is the binding energy and Qq 
utt/c is the Fermi momentum. 

In the strong coupling limit (Qo 
chemical potential tends to zero as 

2 / 3 / h \ 2 / 3 
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FIG. 1: Phase diagram in the plane /i/es, /i/es. The par- 
tially polarised phase (PP), fully polarised phase (FP) and 
vacuum (V) are displayed. The dashed lines correspond to 
the asymptotes discussed in the text; the vertical arrow is 
discussed in the frame of the Local Density Approximation. 



while for the weak coupling limit (Qo ~~ * 00 ^ h 3* es) 
the chemical potential diverges according to 



JL 



h 



(9) 



which recovers the result obtained using a mean-field ap- 
proach. Since the excitation spectrum is gapless, the sat- 
uration line crosses the point h = 0, /a = 0. The vacuum 
line, for which n = 0, corresponds to fi v = —h. The 
phase diagram is shown in Fig. [1] 

We will now calculate the magnetic susceptibility close 
to saturation. Using the equality 



which only diverges in the limit of the empty band 
(Qo = 0) and in the limit of infinite interaction, and 
remains finite otherwise. We want to stress that the cal- 
culations have been performed for the case with fixed 
number of particles, for which the van Hove singularity 
is not present, therefore we do not expect a divergency 
of the magnetic susceptibility at saturation. 



IV. TRAPPED DENSITY PROFILES 

Let us now add a harmonic confinement along the axial 
direction V ex t(z) = muilz 2 /2. If the trap is sufficiently 
shallow then the size of the cloud can be considered much 
bigger than a z and one can treat the system as being lo- 
cally uniform. This allows us to use the Local Density 
Approximation (LDA) in order to write the local equilib- 
rium conditions 



/i(n T (z),n x (z)) = h°, 
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2 2 

mu) v z , 



(14) 



where 



h(n^z), ni (z)) = M"T(*).ni(*))-Mn T (*).ni(*)) ( 15 ) 

are the chemical potentials calculated for the homoge- 
neous system, and /x° = (/i°+/i°)/2 and h° = (£t? — A t °)/2 
are constants calculated imposing the normalisation con- 
ditions for the number of particles in each component, 
i.e. for the total number of particles N and the total 
magnetisation S 



f(\)a(X)d\ = /(A* 



a(X)d\, -B < A* < B, (10) 



-B J -B 

we can rewrite Q (5J| and the energy (J3J) as: 
„ „ n — s 
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dp . (12) 



After solving Q in (|11|) up to second order in the small 
parameter (n — s), we can calculate the energy near sat- 
uration using (|12p. Taking the first derivative of the en- 
ergy with respect to the magnetisation and the number 
of particles at n = s we recover the result for the mag- 
netic field and the chemical potential at saturation 0. 
Taking the second derivative of the energy with respect 
to the magnetisation we obtain the value of the magnetic 
susceptibility at saturation: 
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Note that in our situation the two species are confined in 
the same potential (since we can consider them as being 
two different hyperfine states of the same atom) and so 
the local magnetic field h(z) is kept constant along the 
trap, while the local chemical potential fi(z) decreases as 
we approach the edges of the cloud (see arrow in Fig. [TJ. 
In order to calculate the density profile, we solve numer- 
ically Eq. (|14|) fixing the total number of particles and 
the magnetisation: 
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,- - r n 
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(17) 



R 



i.e. fixing the values for \i° and h°. Here z — zaio / V2a 2 , 
n = n/c and s = s/c. The density profiles for different 
values of the polarisation P = (iVj — Ni)/N are shown 
in Fig. 1 
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FIG. 2: Density profile for polarisation P — 0.7 (top) and 
P = 0.15 (bottom) for N{a\ D /al) = 1. The solid line is the 
total density n(z)/(yN/a x ), the dotted line is the density 
of up spins n-f(z)/(\i r N/a z ), the dashed line is the density 
of down spins ni{z)/{s/N /a z ) and the dashed-dotted line is 
the magnetisation m(z)/(\/N /a z ). We can see the two shells: 
the inner one contains an imbalanced mixture of up and down 
spins, while the external shell is fully polarised. 



We can understand the density profile as follows: when 
moving along the trap, the chemical potential difference 
h{z) remains constant while the chemical potential n(z) 
decreases parabolically and so does the one-dimensional 
local density. The interaction strength j(z) = c/n(z) in- 
creases while we move away from the center of the cloud, 
since it is inversely proportional to the one-dimensional 
local density. The system is gapless, therefore for any 
applied magnetic field h it is in the mixed phase, that is 
to say in a mixture of up and down spins, s ^ 0. That 
explains why moving from the center of the trap towards 
the edges of the cloud we find two different phases: at the 
center we have an imbalanced mixed phase for any mag- 
netic field h, while at the edges, where the two species 
become more repulsive, we find a fully polarised phase. 

In three dimensions the density profiles at weak inter- 
action strength present a very similar qualitative struc- 
ture [23I ]. However, when increasing the interaction 
strength in three dimensions a symmetry breaking oc- 
curs, driven by the competition between the repulsive 
interaction energy and the kinetic energy, and a phase 
separation occurs: the minority component is pushed 



0.2 




0.1 



0.2 

h/(e B N(a 1D /a z C 



0.3 

2, 



0.4 



0.5 



FIG. 3: Top: Radii of the inner (inferior line) and the outer 
(superior line) shells. Bottom: phase diagram in the plane 
P vs. h/e B - For N(a\ D /a 2 z ) = 00 (solid line), N(a{ D /a 2 z ) = 
10 (dashed-dotted line), N(a\ D /a1) = 1 (dotted line) and 
N{a\ D /al) = 0.1 (dashed line). 



to the edges of the trap while the majority compo- 
nent accommodates in the center 



231 ] . In one dimen- 



sion we can calculate the limit of infinite inter-species 
interaction exactly by means of the exact mapping to 
a fcrmionic Tonks- Girardeau gas (two-component Fermi 
gas with hard-core inter-species interaction) for which 
the density profile is known to be equivalent to that of 
N = + N\ bosons with Tonks-Girardeau point- like 
interaction [24( . There is though no phase-separation ex- 
pected in this limit in the one-dimensional configuration, 
which is consistent with our numerical calculations. 



We now calculate the radii of the inner and outer 
shells in the density profiles. The first is given by 
{a 1D /a 2 z )R ln = y/2 {ji°/e 
(ai D /a 2 z )R out = ^2 ([J.°/e B - Vv/cb) 



fj-s/es) and the last by 
We compute 



them as a function of the polarisation for different val- 
ues of the interaction strength (Fig. [3j. We see perfect 
agreement between the numerical results and the limit- 
ing analytical result for free fermions, which is displayed 
in Fig. [3] with solid lines. We also see agreement with 
the limiting case of infinite repulsion: this limit is sin- 
gular, since the interaction between the different species 
mimics the Pauli principle and we can no longer differ- 
entiate between up and down spin particles. In this case, 
the radius of the cloud is y2N a z and is independent of 
the polarisation [24|]. We also display in Fig. [3] the varia- 
tion of the polarisation with the magnetic field for a fixed 
number of particles. The onset of magnetisation occurs 
at h = 0, proving the absence of a gap in the system. 
The slope of the magnetisation at the onset of polarisa- 
tion is sensitive to the rate at which the gas polarises; 
the linear behaviour is a sign of a finite susceptibility 
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[25j . The slope remains close to a constant in the regime 
< h < h s , showing that the polarisability of the system 
is very weakly affected by the extra imbalanced electrons. 
Even though there is a discontinuity in the slope of the 
polarisation versus the magnetic field at the saturation 
point, we do not expect a divergency of the magnetic 
susceptibility at saturation (|f 3[) . 



V. CONCLUSIONS 

In conclusion, we have given an exact solution for the 
phase diagram of the ground state of a one-dimensional 
two-component Fermi gas with inter-component repul- 
sive interactions, for any value of the coupling strength. 
Three phases are localised: a balanced phase that only 
exists at zero magnetic field, a phase with an imbalanced 
mixture of the two components, and a fully polarised 
phase containing only the majority component. We show 
that inside a harmonic trap the cloud presents a double- 
shell structure, with an imbalanced mixture of the two 
components in the center and a fully polarised edge. The 
radii of the inner and outer shells are calculated for differ- 
ent values of the polarisation and the coupling strength. 
We finally calculate the dependence of the magnetisation 
on the polarisation for a fixed number of particles, and 
we show that the susceptibility never diverges. 

This model is experimentally accessible. One can 
cool K atoms to quantum degeneracy, then using 
radio-frequency pulses prepare a spin mixture of the 
atoms in the different spin states \F = 9/2, mp = —9/2) 
and \F — 9/2, mp — —7/2) with different concentrations. 
Adiabatically superposing a 2D optical lattice leads to 
arrays of nearly identically parallel ID traps. 
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